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We calculate the fluctuation spectrum of the shape of a lipid vesicle or cell exposed to a non- 
thermal source of noise. In particular we take into account constraints on the membrane area and the 
volume of fluid that it encapsulates when obtaining expressions for the dependency of the membrane 
tension on the noise. We then investigate three possible origins of the non-thermal noise taken from 
the literature: A direct force, which models an external medium pushing on the membrane. A 
curvature force, which models a fluctuating spontaneous curvature, and a permeation force coming 
from an active transport of fluid through the membrane. For the direct force and curvature force 
cases, we compare our results to existing experiments on active membranes. 
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I. INTRODUCTION 

Lipid membranes are an essential part of any biological 
cell, for example as the interface between the cell and its 
surrounding environment. The cell is actively engaged in 
interactions with its environment which means that the 
membrane is kept out of thermal equilibrium. One of the 
interesting properties of the membrane are its shape fluc- 
tuations, which can influence for instance membrane ad- 
hesion [iHS]. The shape fluctuations must be influenced 
by the activity since it is for example known that the fluc- 
tuations of red blood cells depend on the concentration 
of ATP IT . Artificial vesicles are a valuable model for liv- 
ing membranes as well as an interesting physical object 
of its own |7j. The simplest fluid vesicles are in thermal 
equilibrium. Their fluctuations are purely thermal and 
characterized by the well known Helfrich effective free en- 
ergy [S]. It contains the bending rigidity, the Gaussian 
bending rigidity and the tension. The bending rigidity 
and the Gaussian bending rigidity are material param- 
eters, which depend mainly on the membrane composi- 
tion, while the tension can be interpreted as a Lagrange 
multiplier for a constraint on the membrane area [5j . The 
tension is therefore related to the excess area stored in 
the fluctuations of the membrane shape. Increasing the 
excess area often leads to a response of decreasing tension 
and vice versa. A number of theoretical and experimental 
works on active membranes consider the active contribu- 
tion in the membrane dynamical equations disregarding 
effects on the interplay between tension and excess area 
of the membrane [TDHT3] . The correct physical treat- 
ment of the renormalization of the tension is an impor- 
tant problem since, for instance, the membrane tension 
has been shown to affect membrane endocytosis [13] or 
mcchanosensitive channels jl6H18j as well as influencing 
the motility of some cells jT9]. In this paper we use a 
simple model for an active force on the membrane to ex- 
tend the idea of area conservation to active membranes. 



To do so we compare the passive and active fluctuation 
spectrums of the same vesicle and we calculate the ten- 
sion selfconsistently through a constraint on the excess 
area. We will consider three different physical scenarios 
for the force that have been studied in the literature: di- 
rect force [m [20] , curvature force piTHT^ and a case we 
will call permeation force [101 EI] • For the direct force 
we will consider the example of adenosine triphosphate 
(ATP) activated cytoskeleton of red blood cells, [22H21] 
and references therein. For the curvature force we will 
comment on experiments with giant unilamelar vesicles 
having active ions pumps included in their membranes 
PUI [551 US] ■ We will give order of magnitude estimates 
of the effects for the permeation force case. 

The outline of the paper is as follows: In Sectionjlll we 
recall the dynamical equations governing the fluctuations 
and the area conservation for passive vesicles. In Section 
III we add an active contribution to the dynamic equa- 
tions for the vesicle and give the corresponding changes in 
the general expressions for the fluctuations and the area 



conservation equations. In section IV we discuss the di 



rect force case; in Section [V] we introduce the curvature 
force case, and in Section ]VI| the permeation force case. 
Finally, we give our conclusions and outlook in Section 

ivni 



II. EQUILIBRIUM FLUCTUATIONS 

We begin by recalling the equations of motion for the 
passive fluctuations of a quasi spherical vesicle taking into 
account the area and volume constraints |9j . We take the 
following effective free energy for our membrane 



F = j dA {2kH^ 



(1) 



The first term in this equation is the classical Helfrich 
energy with k being the bending rigidity, H the local 
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mean curvature and the integration being over the whole 
area of the membrane. The second term is a Lagrange 
multipher for the area: the tension S wiU be determined 
selfconsistently such that the total area is conserved. For 
the quasi-spherical vesicle we will write the radial coor- 
dinate of the surface as 



R{e,^,t)=Ro{l + u{0,(j),t)), 



(2) 



where k = K/{kBT) is a reduced bending rigidity. We 
now calculate the tension ath as in [9] . We start by writ- 
ing down the area and volume of the vesicle as a function 

of the Uim's 



A - i?2 4^ 1 



^00 



where Rq is the radius defined through the fixed volume 
V of the vesicle hy V = 47r_RQ/3, and 9 and 4> are the 
polar and azimuthal angle. We expand the u{9,(f),t) in 
spherical harmonics as 



(3) 



1=0 m=-l 



where 1^ is a cut off value for /. The equation of mo- 
tion for the uim (t) is the normal component of the force 
balance equation of the membrane. Detailed derivations 
have been given elsewhere, [11] or [27] for example. We 
only state the result here 



^^/™(^) = -j^Eiia,^)uUt) + ^rit), (4) 



where 



r, = 



4^3 + 6P~1' 



Ei{a) = {l + 2){l-l)[l{l + l) + a] 



(5) 



(6) 



The left hand side of Eq. Q is the friction of the sur- 
rounding bulk fluid on the membrane with 77 being the 
viscosity of the fluid and the dot denoting a time deriva- 
tive. The first term on the right hand side is the elastic 
restoring force coming from the Helfrich free energy Eq. 
^ where ath = Si?^/ K is a dimensionless tension. The 
last term is the thermal noise giving rise to the thermal 
fiuctuations. The properties of this thermal noise are 



= 



(7a) 



'S{t-t'), (7b) 



where the star denotes complex conjugation, fc^ is the 
Boltzmann constant, T the temperature of the surround- 
ing fluid. Sab the Kronecker delta and S{t) the Dirac delta 
function. Here and throughout the paper the brackets 
mean that we average over noise realizations. From ei- 
ther these equations or the free energy we can calculate 
the u;,„ fluctuations to be 



J2 E \ni,nf (1 + 1(1 + l)/2) + 0iufJ 

; = 1 m=-l / 



(9) 



/47r, 

£ E l"^™l' + o(uL) • (10) 



Having a fixed volume V = 47ri?o/3 fixes the amplitude 
Moo a-s a function of the other uim'- 



uoo 
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(11) 



The excess area, which we will define as A = (A 
4nR'^)/Rl, is then 



In-, I 

^-E E \ui,n\'^(l + 2)(l-l). 



(12) 



1=2 \ 



Note that we did not include a factor of Att in the def- 
inition of the excess area for convenience and that the 
I — 1 mode, which corresponds to translation, cancel. To 
obtain an equation that gives the tension as a function of 
the excess area we will consider a thermodynamic limit 



where we can replace |ii;„J^ by ( |u;„iP) in Eq. ([12 1 

' \ / th ' — ' 

It then reads 



in 



21 + 1 



1=2 



1(1 + 1)+ ath 



(13) 



This gives a self consistent equation for the passive ten- 
sion ath that has been studied in 



III. FLUCTUATIONS IN THE ACTIVE CASE 

To model the activity we add an active random con- 
tribution in the force balance equation Eq. Q. The 
equation of motion in the active case is then 



^Ui,n(t) 

1 1 



Ei(aMn.(t)+^l:^\t) + ^rj(t), (14) 



(tKm(t)) 



1 



th 



iiEi(ath) 



(8) 



where we introduced a new subscript in the symbol for 
the tension Ca in the active case, indicating that it can 
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be different from the tension ath in the passive case, and for the area conservation in both active case, Eq. (191, 
^im (0 the noise coming from the active processes. We 
wiU take this noise to have the foUowing properties 



and passive case, Eq. ( 13 ), we can see that da > Cth with 



' exp 



(15a) 



(15b) 



where Ta is a characteristic correlation time of the ac- 
tive process studied. It could be on the order of the 
average time for a pumping cycle for an ion pump for 
example. Xa gives the strength of the noise. It has the 
dimension of a force per unit area squared (N^m~^). xi 
is a dimensionlcss quantity that carries the wavenumber 
dependency of the noise. By spherical symmetry it can 
only be a function of I and not to. From the Langevin 
equation ( 14 1 and the noise properties Eqs. ( 15 1 and Q 
the uim fluctuations can be calculated as 



\Ulm\^) = (l+Xa ^^'_ , (16) 



K^i(c^a) 



where 



Xa = X; 



a 7 rri 5 
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(17) 



(18) 



Xa is a non-dimensional measure for the strength of the 
noise and tm is the correlation time for mode I of the 
shape in thermal equilibrium. The fluctuation spectrum 
is similar to the passive fluctuation spectrum with the 
active tension CTa replacing the passive tension dth and 
an additional positive contribution proportional to the 
strength of the noise Xa- The effective equation of area 



conservation (12 1 (with \ui 



2»A = Ei^jTI)T7 




Equation ( 19 ) is a central result of this paper. It relates 



the excess area A and the tension to the parameters 
describing the activity: Xa, xi and t^- We see that the ef- 
fect of the activity is twofold: it give a direct contribution 
to the fluctuation spectrum Eq. ( 16 1 and it indirectly 



renormalises the tension through the area constraint Eq. 
pl. 



Tension renormalization 



Using Eq. ( 19 1 we can deduce some general properties 
of the tension in our model. Comparing the equations 



CTa = CTth only if Xa = 0. Because we have fixed the area of 
the membrane then the additional noise due to the active 
forces will always lead to a counterbalancing increase in 
the tension. Also note that in our model the total tension 



that appear in the fluctuation spectrum, Eq. ( 16 1, is fixed 
selfconsistently entirely through Eq. (19). This implies 



that any direct contribution to the tension in the force 
balance equation ( [T4| ) has no effect on the fiuctuation 
spectrum, since it will be canceled by a counterbalancing 
tension contribution to keep the excess area conserved. 
In particular the direct contribution to the tension from 
the dipole moment of a distribution of force due to an 
active protein found in [2H! does not affect the fiuctua- 
tion spectrum in the present model of an incompressible 
membrane (see [29] for the case of a slightly elastic mem- 
brane). However, the quadrupole contribution found in 
P5] will act like the curvature force discussed here in Sec- 
tion [V] and thereby it indirectly renormalizes the tension 
through the area constraint. 

We now investigate the relation between the active ten- 
sion (Ta, the excess area A and the other parameters of the 
system using Eq. ( [l9| in different limiting cases. First 
we distinguish two limits: long and short active correla- 
tion time. For long correlation time, formally Ta S> tm , 
Eq. ( 19 1 becomes 



A = ( 1 

6 -f CTa 



E 

J=3 



X2 
"E2{<Ja 

21 + 1 



l{l + l)+aa 



xi 



Ei{aa) 



where A is the reduced excess area 



A = 2Ak = 2A, 



For a short correlation time, formally <C t 



, (20) 

(21) 
we get 



A 



1 + XllTX2 



6 



E 

1=3 



55, 
2/ + 1 



CTa + Z(/ + l) 



where f reads 



{l + XaTXiTi), (22) 



(23) 



Note that this case corresponds to a timewise Dirac delta 
correlated active noise, i.e., a white noise. It can be ob- 
tained as f — > while Xa ~> oo such that x^t is finite. 
In both cases, Eq. (20 1 and Eq. (22), we have isolated 
the I — 2 term as this term diverges when CTa — >■ ~6. In 
the last term of both equations we will approximate the 
sum by an integral while the first term will be kept as it 
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is. Using this approach we can invert Eqs. ( |20[ ) and (22) 
in order to obtain a relation between (7a and the reduced 
excess area A. Similarly to the treatment of the passive 
case in [S] we also break up our results depending on the 
value of (Ta: cr& -6, 1 < cia < and 1 < ^ 
After expanding in the appropriate limit of Ca and Z,,! we 
keep only the dominating term proportional to Xa as well 
as the dominating term that does not depend directly on 
Xa in Eqs. (20) and (22 1. This procedure gives the re- 
sults presented in Table l] for each of the three cases of 
force that we will motivate and consider in the following 
sections. These results are the essential findings of this 
paper. They give the behavior of the tension in terms of 
the active parameters Xa and Ta in the different physical 
cases we consider here. These cases, classified by their de- 
pendency of xi on the wavenumber /, are: a;; ~ 1 (direct 
force), xi ~ P (permeation force) and xi ~ l"^ (curvature 
force). 



IV. DIRECT FORCE 

The case Xi = 1 corresponds to a noisy force that is 
being applied locally on the membrane, i.e., a direct ran- 
dom force contribution in Eq. (14 1. This kind of noise 



have been studied previously, mainly to model an active 
coupling between the cytoskeleton and the membrane in 
red blood cells [211 [Ml [SI]- In this case it represents 
the punctual force applied on to the membrane (or more 
precisely, a release of tension) when ATP is used in the 
process of cutting a link between an actin filament and 
the membrane. The amplitude of the active correlation 
can be estimated as 



Xa 



El 
Ro 



(24) 



where Fq is on the order of the force applied on the mem- 
brane by one active center and p is the concentration of 
active centers. In the cytoskeleton case p will be related 
to the average density of links activated by ATP and 
hence it will depend on the ATP concentration [53] . The 
active correlation time Ta should be on the order of the re- 
lease time of an actin filament. Note that the direct force 
case requires an external source of momentum for New- 
ton's third law to be obeyed; the case cannot be applied 
to proteins included in the membrane without external 
attachment, for example ion pumps. 



[T] we have shown the renormalisation of the tension as a 
function of the strength of the activity with some typical 
parameter values for a red blood cell. Here we compare 
this result to another model [551 01] • This model intro- 
duced a confining potential that act on the membrane in 
order to model the attachment of the cytoskeleton to the 
membrane. In addition, the calculation was performed 
for an almost flat membrane with infinite frame area. 
In the spherical limit we consider in this paper we can 
get the flat case limit by keeping only the highest power 
of I and then replacing I by qRo where q is the Fourier 
transform variable. In this limit Ei, including a uniform 
potential, becomes 



kEi 4 2 

— ^ -> + erg +7. 

Kr, 



(25) 



The parameter 7 is the strength of the conflning poten- 
tial. There were no constraints on the area but a effective 
(increased) tension was considered to come from the con- 
finement of the membrane |22j . We can give the following 
two remarks here comparing to our approach. First, in 
our model a laterally uniform harmonic confining poten- 
tial over the membrane will decrease the tension com- 
pared to a case where there is no confining potential (no 
attachment to the membrane) because the uniform po- 
tential tends to decrease the low wavelength fluctuations; 
hence the absolute value of the tension will be diminished 
by the confinement. Secondly, the increase of the tension 
in our model is due to the effect of the noisy activity 
through the area constraint. Note that with increasing 
ATP concentration the strength of the activity Xa will 
increase while the confinement will decrease because the 
cytoskeleton will have less attachments to the membrane. 
Thus taking into account the area constraint leads us to 
conclude that, in our model, when the ATP induced ac- 
tivity increases the tension increases. This contrasts the 
modelling of [551 US] , where a non-homogenous potential 
due to the presence of the cytoskeleton induces a positive 
contribution to the tension, which implies that the ten- 
sion should decrease with increasing ATP concentration. 
We will show in the next section that our result of an 
activity induced increase of the tension is not in contra- 
diction with the decrease of the fluctuations observed in 
ATP depleted red blood cells. Also note that the tension 
we calculated here is only partially related to the force 
needed to pull tethers out of cells, as the attachment of 
the membrane to the cytoskeleton directly contribute to 
this force as discussed in 1551. 



A. Tension renormalisation 

The activity will modify the tension through the area 
constraint as described in Table |T1 We begin our discus- 
sion of the application to the cytoskeleton by pointing out 
that the ATP induced activity, by attempting to increase 
the membrane fluctuations, increases the membrane ten- 
sion through the area constraint. Hence the tension in- 
creases when the ATP concentration increases. In Fig. 



B. Fluctuation renormalisation 

In Fig. [2] we compare the passive fluctuation spec- 
trum, Xa — 0, with the corresponding active one for 
Fq ~ 0.2 pN and the same parameters as in Fig. [ij We 
can see that the activity increases the large wavelength 
(small I) fluctuations while decreasing the small wave- 
length (large /) fluctuations; following the general argu- 
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TABLE I: The different behaviors of the tension obtained for the direct force, curvature force and permeation force in the hmit 
of short and long correlation times for the active noise. 



a^ = -6 + \ 1 
2A 



aa = -6+^(l + 4xa?r2) 



aa = -6 + ^ (1 + Xaf Fa) 



5 A 55 _ _ 



(Ta = exp 



A + Xa/4' 



1 + Xa/4 



(Ja = exp -A + 



4 6 



A = ln('^U^flnr'^ 



A = In 



4 A/Oa 



A = ln( +16^ 

. (Ja / fa 



(Ta = exp -A + 8XaT/ 



1 < li < (Ta 



.a^i| l + ^l+iAx. 



^ + 40 



1+ /l + 4%Aln('^ 
2A \ V li \ 4 



^ _ A , Xa^ 



O-a = ^ I 1 
2A 



1 + 64A^ 



O-a = ^ I 1 + gXa"n 



ment laid out in the supplementary material we can show 
that, in the direct force case for dth > 0, there exists a 
crossover I = Ic such that the activity always increases 
the fluctuations for I < Ic and decreases the fluctuations 
for I > Ic- In the inset of Fig. [2] we have plotted the 
amplitude of the fluctuations in real space 
which can be evaluated as 



{\u{d,cl>)\')^{\uf)^]-f2 J2 (I 



Ulr. 



(26) 



1=0 m=~-l 



This amplitude is independent of 9 and (f> due to rota- 
tional symmetry. From the inset we can see that the real 
space fluctuation amplitudes can increase by a factor of 4 
due to the activity even though the tension is increased 
and the excess area is conserved. Such an increase of 
the fluctuations has been observed in [33j by varying the 
ATP concentration in red blood cells. The effect of the 
direct force activity is to increase the large wavelength 
fluctuations (even though the tension is increased) while 
decreasing the fluctuations in an intermediate range of 
wavelengths. This increases the observed real space am- 
plitude of the fluctuations while still keeping the total 
excess area unchanged. 



V. CURVATURE FORCE 

The curvature case has also been studied previously. 
It is presented as a random active contribution to the 
spontaneous curvature in [30'. It can also be seen as 
a random contribution to the quadrupole moment of a 



microscopic force density due to active transmembrane 
proteins [TOl [28] . Note that any second order derivative 
of a field on the membrane surface, in the force expression 
in real space, would give a "curvature force" noise in our 
linear theory. Such a noise will lead to the scaling xi ^ 
at large I. For the subdominant behavior we will follow 
the behavior found for the quadrupole contribution in 
[TT] and therefore choose xi = {I + 2)^(1 — 1)^/4 for this 
case and estimate the strength of the noise by 



Xa 



P 



F2 

i?3 



(27) 



where F2 has the unit of force times length squared, p 
is the concentration of active proteins per area of the 
membrane and Ta is a time scale on the order of the 
pumping cycle time for an ion pump for example. 



A. Micropipette experiments 

In [lOl [25] vesicles containing active ion pumps have 
been studied experimentally using a micropipette aspira- 
tion technique. In these experiments a vesicle is sucked 
partially into a micropipette while measuring the change 
in excess area and the tension applied on the membrane. 
In both cases they activated the pump externally (by 
shinning light on them or by injecting ATP). They com- 
pared the relation between the excess area and the ap- 
plied tension for the same vesicle with and without activ- 
ity. For a passive vesicle the relation between the tension 



and the excess area is given by Eq. (13). In this case 
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15 
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arithm of the tension still holds albeit with a decreased 
prefactor. This decreased prefactor was assigned to an 
increased effective temperature T^g. Using the equations 
in Table U we can relate the excess area and the tension 
in the curvature case for Ta » t^^ and 1 ^ CTa ^ ZH, as 
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FIG. 1: The relative renormalisation of the tension as a func- 
tion of the strength of the activity in the direct force case. The 
horizontal axis represents Xa defined by Eq. ( |24[ | where Fq 
go from to 0.2 pN. The big circles represent era's calculated 
from Eq. ( 19 1 while the dotted line is (Ta calculated from the 
formula in Table \U in the direct force case for 1 ^ (Ta ^ 

(2) — 

and Ta 3> tm . The paramaters used are: Rq ~ 5 ^im, 
77 ~ 9 X lO'^kg/ms, ksT ~ 4 x 10"^M, k ~ AksT, 
Ta — 0.1 s, assuming a typical distance between two active 



centers of 100 nm we have p ~ 10 m" 



^, and taking a mem- 
~ Ro/d ~ 1000. We 



brane thickness of ci ~ 5 nm we get k 
also chose A/47r ~ 3% from which we found ath — 529 by 
numerically solving Eq. (13 1. 



0.001 - 



10" 



A 



10" 



s 

V 10" 
10"' 
10"' 



15 : 
V 10 

/ s' 



0.00.5 1.01.5 2.02.5 3.03.5 
10' i„ 




500 1000 



FIG. 2: Two fluctuation spectrums for the direct force case of 
the same vesicle with no activity (xa = 0, continuous line) and 
with Fo — 0.2 pN (dashed line). Inset: the relative increase of 
fiuctuations in real space as a function of the active amplitude. 
The parameters used are the same as in Fig. [l] 



a small increment of the excess area dA is related to an 
increment of the logarithm of the tension d(lncrth) in the 
case 1 < cTth < by 



-dA 



-d(ln(Tth) • 



(28) 



For the active membrane they found that this linear re- 
lationship between the change in excess area and the log- 



2— dA 



Xa 



d(ln(Ta) • 



(29) 



Thus in our model an effective temperature can be de- 
fined as Tcs/T = 1 + Xa/4. The difference between 
our model presented here and the one proposed in [TD] 
is that we look at a noisy force distribution with zero 
mean quadrupole moment while [10] had a non-noisy 
quadrupole force for each protein, which influences the 
fiuctuations through a coupling to the protein density 
and the curvature of the membrane. Both approaches 
lead to the linear relationship of Eq. ( 29 ) , with the cor- 



respondance that F2 — 2wVa where w is defined as the 
membrane thickness in |TOj and Va is labeled the force 
dipole. In the case of bacteriorhodopsin they found ex- 
perimentally that Teff/T ~ 2. By taking k ~ lOfcsT 
and p ~ 10^^ m~^ we find from this the quadrupole mo- 
ment F2 ~ 2.6 X 10~^^ Jm. It should be noted that the 
effective temperature T^g defined here for the curvature 
case is related to the increase of the fiuctuations at small 
wavelength of the fiuctuation spectrum. Indeed taking 
the limit of the wavenumber Z — > 00 in Eq. ( 16 1 we get 



{\uir. 



ksT 



(1 



4'^ kI^ ' 



(30) 



showing that the same effective temperature appears in 
the small wavelength region of the fiuctuation spectrum. 



B. Videomicroscopy experiment of fluctuating 
vesicle 

Another experiment on artificial active membranes was 
performed in [26j . Here both the active and passive fiuc- 
tuation spectrum of a giant unilamellar vesicle contain- 
ing the proton pump bacteriorhodopsin was recorded by 
videomicroscopy, without the vesicle tension being con- 
strained externally. Before discussing their results we will 
lay out what we would expect in terms of the present the- 
ory: The short wavelength effective temperature found 
in the micropipette experiment should still be present in 
this free vesicle experiment. Thus the active contribution 
should still increase the short wavelength (large I) fluctu- 
ations compared to the passive ones. For the membrane 
to conserve its excess area it needs to compensate for the 
increase of the small wavelength fluctuations through an 
increase of the tension such that the large wavelength 
fiuctuations are decreased. This can be seen in Fig. [3] 
where the expected active spectrum (dashed line) is com- 
pared to the passive one (full line). In fact for the cur- 
vature case it can be proven that the active fluctuation 
spectrum is always larger than the passive one for small 
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I and then become smaller at values of I above a critical 
I = Ic (see the supplementary material for a calculation 
of the IcS for small Xa)- 



If we take the limit of long correlation time: Ta 
of Eq. ( 16 ) for the curvature force case we get 



> t 



(2) 



\Ulr. 



kEi (o-a 



Xa {l + 2)^{l-ir 

4 Ei{<7,) 



(31) 



The flat limit of this equation, i.e., ignoring subdomi- 
nant terms in I e.g. {I + 2){l — 1) ^ P , is the equation 
used to fit the experimentally obtained fluctuation spec- 
trum in if we evaluate Xa using Eq. (27). This is 



true even though we have considered a shot noise with 
zero mean while they have considered a mean contri- 
bution coupled to the proteins density difference (sim- 
ilarly to the discussion in the last subsection). Or in 
other words: a very long correlation time for the noise is 
equivalent to a constant mean contribution in this case. 
Experimentally they found a significant increase of the 
fluctuation spectrum at low I values when the proteins 
were activated, which they interpreted as an active re- 
duction of the tension from Sth — 4 x 10^^ Nm^^ in 
the passive case to Ea — 5.3 x 10^^ Nm^^ in the active 
case. This gives Cth = 940 and Ua ~ 130 for our unitless 
tensions. The corresponding active fluctuation spectrum 
calculated from Eq. ( [3T| ), with the value Xa = 4 taken 
from the fit to the micropipette experiment and tia = 130, 
is plotted in Fig. [3] as a dot dashed line. 



C. Discussion 

The lowering of the tension is necessary to explain the 
increase of the large wavelength fluctuations observed in 
|26j and the increase in effective temperature is necessary 
to explain the micropipette experiment [TU]. But this 
gives an active fluctuation spectrum consistently larger 
than the passive one for all wavenumbers. These results 
suggest that the excess area is not conserved in these 
experiments. We can suggest two different mechanisms 
in order to explain this increase: The first explanation 
would be that the increase of the excess area is due to 
the elastic properties of the lipid bilayer itself and a very 
large negative contribution to the tension in the active 
case. In this case it is the lipidic part of the membranes 
area which is expanded due to the additional negative 
active contribution to the tension. However, the active 
contribution to the tension need to be very large in mag- 
nitude, at least on the order of 10^'^Nm^^, to give a 
visible effect on the fluctuation spectrum due to the high 
compressional modules of a lipid membrane, see |29,. 

Another explanation could be, for these experiments 
using bacteriorhodopsin, that the area the proteins take 
up in the membrane increases when they are active (at 
least on average in time), since this protein is known to 
periodically expand during its active cycle [34j . In this 
case the area per lipid head in the lipid region of the 
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FIG. 3: Three fluctuation spectrums for the curvature case. 
The full line corresponds to the passive case, Xa = 0, taking 
the passive tension obtained in 26 : Sth — 4 x 10~^N m"'^ 
leading to ath — 940. The dashed line corresponds to the 
active case with Xa = 4 and a tension = 3.8 x 10* (giving a 
conserved excess area) . The dot dashed line is the fluctuation 
spectrum using values mimicking the experiments, Xa = 4 and 
CTa — 130, see the discussion in the text. The vertical dotted 
line denotes the approximate maximum wavevector observed 
in [26]. We used the parameters: k ~ WksT, p ~ 10^^ m~^, 
i?o ^ 10 /im and l,n = Ro/d ~ 2000. 



membrane would be unchanged and the increase of the 
excess area would be due to the expansion of the proteins. 
This would allow the membrane tension to decrease. 

Assuming Xa — 4 we can evaluate the excess area Aa 
required to obtain Ca = 130. Using crth = 940, 1^, = 2000 
and K — lOfc^T we find for the passive case Ath/(47r) ~ 
0.03. From the result for the tension in Table ll] for the 

('2') 

curvature case with Ta 3> tm and 1 ^ (Ta ^ obtain 



1 



Xa 
4 



log 



gth 
CTa, 



Xa 

T 



(32) 



With (Ta = 130 and Xa - 4 we find Aa/(47r) ~ 0.07. One 
can evaluate the corresponding increase of the area per 
proteins, a^, between the passive and active case. For 
a protein concentration of p ~ 10^^ m~^ we find ~ 
4 X 10^^^ m^. This correspond to a ^ 10% increase of 
the active radius, if wc take a protein to occupy a disk of 
radius Vp ~ 2.5 nm in the passive case. 



VI. PERMEATION FORCE 

We will show here that a contribution of the form 
xi ~ P can be justified through a model of the activ- 
ity with transfer of fluid between the two sides of the 
membrane. In Eq. ( 14 ) we used a no slip boundary con- 



dition to obtain the hydrodynamic stresses on the mem- 
brane from the surrounding fluid. The radial component 
of this boundary condition reads 



R{e,ci),t)-vU^{e,c^,R,t) = o, 



(33) 
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FIG. 4: The relative increase in tension as a function of the 
mean volume transfered through the membrane per unit time 
for the permeation case. The dotted line is calculated using 
the equation for Ta S> fm' in Table |l] while the big circles 
are numerical calculations directly from Eq. (19 1. Inset: the 
first 30 modes of the fluctuation spectrum. The blue line 
is for the passive case while the red dashed line is for — 
1.5 X 10~^®m^/s. We used the parameters k ~ lOfcsT, p ~ 
lO^'^m-^ Ro ~ 10 pm, 1^ = Ro/d ~ 2000, ~ 10"^ kg/m s 
and (Tth ^ 940. 



where R is the radial position of the membrane and Wflui^ 
is the radial component of the fluid velocity evaluated at 
the membrane. To obtain the permeation force case we 
add an active contribution to the boundary condition. It 
then reads 



coefficient Ap and a force Fa was pushing fluid through it. 

Our parameters are related to theirs as pV^ ~ p(ApFa)^ 
where p = p is the average pump density on the mem- 
brane. They neglected this contribution in their effec- 
tive temperature expression for bacteriorhodopsin be- 
cause the water permeation coefhcient through the lipid 
bilayer is very low, Ap < 10~^^ m'^/(N s); however this 
case could give a significant effect for pumps if they could 
provide a high enough flow through the membrane. Us- 
ing Eq. (361 in Eq. ( 16 1 we found the same equation for the 



fluctuation spectrum as 121^ for a zero mean active flow 
in the flat limit. Unfortunately we are unaware of any ex- 
periment relating fluctuations to fluid pumping through 
the membrane. For a detailed study of the fluctuation 
spectrum in this permeation force case we refer to the 
supplementary material. 

As an example we will evaluate the flow, Va, required to 
obtain a signiflcant renormalisation of the tension using 
the parameters of the experiment described in Fig. [3] We 
have plotted the renormalisation of the tension, in the 
limit of long correlation time, as a function of the mean 
volume per second per protein in Fig. [4] With the 
chosen protein concentration we start to have an effect for 
the renormalisation of the tension around 5 x 10~^" m^/s 
which roughly corresponds to 10^ water molecule per ms. 
This order of magnitude for the flow can be obtained 
for the passive channel aquaporin |35j for instance. But 
we believe it is unlikely that active pumps can attain 
that level of flow through the membrane. Note, however, 
that decreasing the tension will increase the effect of the 
activity for a fixed Xa- 



R{e, 0, t) - vi,^,^{e, 0, i?, t) = vi{e, t), (34) 

where is the random active contribution of zero mean 
and with the following correlation after expanding in 
spherical harmonics 



{<,l^{Q)vl,Vm'{t))=p{t 



7?2 



■ exp 



1^1 



]35) 

Here p is the mean density of active centers and Va is 
on the order of the volume transferred through the mem- 
brane per active center per unit time. We have assumed 
that the active centers are operating independently in 
Eq. ( 35 ) . Using Eq. ( 34 1 instead of the no slip boundary 



condition we obtain a dynamical equation of the form 
of Eq. ^ with et'W = Vvl,^{t)l{TiR^). 
parameters Xa and xi then read 



The active 



Xa 



Rl , 



and xi = {Ti) 



(36) 



which gives xi ~ P at large /. 

In |10j a similarly active permeation term was consid- 
ered for the force dipole model using Darcy's law of per- 
meation. The membrane was taken to have a permeation 



VII. CONCLUSION 

In this paper we have investigated the fluctuation spec- 
trum of vesicles with noisy active inclusions. We used a 
simpler model for the activity than previous studies ^iDl - 
I12j . since our noise had vanishing mean, but instead we 
took into account the high compressional modulus of a 
lipid membrane by an area constraint. We investigated 
the effect of the area constraint on three physical models 
for the activity. We showed that our model for the di- 
rect force case agree well with previous experimental and 
theoretical work on red blood cells; albeit with a mod- 
ified interpretation for the increase of the tension. For 
the force curvature case we compared our model to the 
experiments of (TUl HHj and we showed that even though 
our model recover the effective temperature found in the 
micropipette experiment, the data for the free vesicle 
containing active bacteriorhodopsin seems at odds with 
the conservation of the excess area. This suggests that 
bactheriorhodopsin might increase the excess area of the 
vesicle when the proteins are active or that the force cur- 
vature case does not completely describe those experi- 
ments. We also investigated a permeation force case and 
gave an order of magnitude of the flow required to ob- 
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tain a significant tension renormalisation for a specific 
example. We hope that this study wiU motivate further 
theoretical and experimental works on active biomem- 
branes. 



Appendix A: Supplementary material 

In this supplementary material we investigate a large 
part of the phase space of the possible behavior of the 
active fluctuation spectrum, compared to the passive one. 



1. General formalism 

We study the behavior of the active fluctuation spec- 
trum (xa > 0) compared to the passive one (xa. = 0) for 
the same vesicle taking into account the area constraint. 
To do so we study the sign change as a function of I of 



X/l |2\ - d^' 



m /a 



(SI) 



ml /th 



where {\uim\^)a is the active fluctuation spectrum calcu- 
lated in the main text, Eq. (16), and )th is the 
passive fluctuation spectrum Xa = 0, Eq. (8) in the main 
text. After a little algebra 5(|Mzm| ) reads 



S{\uir, 



l{l + 1) + cra 



{(^th - fa) 

£^;(fth) 



\l + 2){l~l) Ei{a,) r, + t 



(0 



(S2) 



The first term in the parenthesis is negative and con- 
stant in / while the last one is positive and depends on /. 
The factor in front of the parenthesis is alwajs positive, 
hence it does not change the sign of (5(|uim| ). For the 
area to be conserved there is at least one sign change of 
S{\uim\'^) in the accessible I range: if the active fluctua- 
tions are greater than the passive fluctuations at some I 
value it must be lower at some other I value to have the 
same excess area and reciprocally. In the rest of this sec- 
tion we lay out our general strategy to investigate these 
cross-overs. We will then specialize to the physical cases 
studied in the main paper in the next three sections. 
The first thing we can do is to check if the terms in 



the parenthesis in ( S2 ) is monotonously increasing or de- 



creasing in I by studying the sign of its derivative given 
the physical conditions (Tth > —6, Cth < (Ja, 2 < / < l^, 
Xa > 0, K > and Ta > 0. If for these conditions the 
derivative does not change sign over 2 < Z < then the 
function in the parenthesis of (S2 ) is monotonous in I. If 



this is the case then we can conclude that (5(|u;m| ) has 
only one sign change in the region 2 < I < i.e. there 
is only one crossover between the active and the passive 
fluctuation spectrums. We can then deduce if the active 



fluctuations are larger/smaller than the passive fluctua- 
tions before and after the crossovers depending on the 
sign of the derivative. 



If the function in the parenthesis of ( S2 1 is not 



monotonous there might be one or more crossovers de- 
pending on the parameter values. In order to study 
these cases we consider two limiting regimes: long and 
short active correlation times. For long correlation time, 

.(2) 



Ta > (S2) reads 



1 



(Cth - fa) 



, ~ Xl Ei{ati,) \ 



while for short correlation time, Ta ^ t[t"' , it reads 
1 



S{\uim\ ) 



1(1 + 1) + 



(dth — Ta) 



{l + 2){l^l) 



Eiiati,)ri\. (S4) 



Then we can again try to differentiate the terms in paren- 
thesis of (S3) and (S4) to check if they are monotonous. 



If this is the case then the above reasoning apply (there 
is only one Ic etc.). 

If the function is not monotonous then we resort to an 
expansion in small Xa, the precise meaning of which de- 
pends on the case studied. We expand the active tension 
CTa as 



fth 



(S5) 



where e <C Cth is the difference between the active and 
passive tension, e is related to Xa. through the excess area 
constraint. To first order in both e and Xa we can rewrite 
Eq. (19) in the main text, in the limit Ta 3> tm , as 



2/ + 1 



[?(/+ 1) +(Tth] 



Xl 



. 21 + 1 

^^'^/(; + l)+'Tthi?((fth)' 



(S6) 



and in the limit Ta <C t 



we get 







1=2 



21 + 1 



[/(/+ 1) +CTth] 



la 



21 + 1 



^ ?(Z+l)+f7th 



We have assumed that the passive tension ath satisfies 
Eq. (13) in the main text with the same excess area as 
in the active case. Xa = implies e = 0. The relation 
between the two then reads 



ClXa 



(S8) 



10 



where ci is the first order coefRcient in an expansion of 
e in Xa <tC CTth/ci- Using the Mac Laurin formula we can 
calculate the sum in (S6) and (S7|. We then consider 
three different cases: dth — >■ —6, 1 ^ ath ^ ^■'^d 
1 ^ ^ (Tth and we expand ci in the appropriate 
power of (Tth and 1^- The lowest relevant order of ci are 
shown in Table [Til 

Next we search for the I — l^s where the magnitude 
of the active fluctuations cross the passive ones. To do 
so we solve the equation i5(|u;m|^) = to first order in e 

(2) 

using (Ta = CTth + CiXa- For tin the Ic are solutions 
of 



and assuming Ua ^ '^^ get 



Cl 



xi 



{i + 2){i-iy 



while for ^ im"'' we have 



Cl = XlTi [l{l + 1) + (Tth] • 



(S9) 



(SIO) 



Note that the crossover values of Ic are independent of Xa. 
and r for small Xa- The Ic values we found are shown in 
Table [TTl We can get the lowest order of ci in an expansion 
in (Tth and /m directly from the equations in Table 1 of the 
main text. Unfortunately the lowest order is not always 
sufficient to determine a non trivial Ic value. 

In the cases where the term in the parenthesis of 
S{\uim\^) was not monotonous and Xa > Cth/ci we have 
calculated the spectrum numerically (in addition to the 
first order expansion). Note that when Xa ~^ c« the 
excess area conservation, Eq. (19) in the main text, 

tells us that the ratio Xa/^a for '''a >• ^m' and Xa/^a 
for Ta ^ im"''' is constant. This can be observed also in 
the last column of Table 1 in the main text. Inserting this 
result in Eq. (16) in the main text we can see that there 
is a point, when Xa is sufficiently large, where the fluctu- 
ation spectrum is not changing significantly anymore as 
Xa increases (note that the tension is still increasing 
with Xa)- We can then investigate the fluctuation spec- 
trum beginning with a small value of Xa ^ fth/ci and 
then increase its value until the fluctuation spectrum do 
not change signiflcantly anymore. 

We now investigate each specific case following the 
method we laid out in this section. 



Curvature force 



In the curvature case the term in parenthesis in ( S2 ) is 
always increasing in I. 6{\uim\ ) goes from being negative 
to being positive as I increases, i.e., the active fluctuation 
spectrum is smaller than the passive one at low / and 
larger at larger I with only one crossover between the 
active and the passive fluctuation spectrum. It can also 
be seen (less generally) this way: if we look at the active 
fluctuations, Eq. (16) in the main text, for the curvature 
case and take the limit I ^ oo, for sufficiently large Ta 



\Ulr. 



:(1 + Xa 



4(^2 + (Ta) 



), (sii) 



where we can see that the activity increases the large I 
fluctuations. The only way the membrane can compen- 
sate for this increase is by making the low / fluctuations 
smaller through an increase in the tension. This is also 
the case in the limit Ta <^ tm"'' albeit that we get a dif- 
ferent I dependency. In the limit Z — >■ (30 we get 



(Knr)a-^^^,(^^(1+Xar{g). 



(S12) 



This increase of the large wavenumber fluctuations can 
be seen as the activity dragging the excess area toward 
large I. In the curvature case the behavior of the ac- 
tive fluctuation spectrum relative to the passive one are 
similar in between short <^ im"' and long Ta 3> 
correlation time. The form of the fluctuation spectrum 
changes, however, as the I dependency of the active part 
of the fluctuation spectrum changes in between short and 
long active correlation time. 

We show an example in Fig. [5] for all the different 
limits of the correlation time and the tension for purposes 
of comparing to the direct and permeation force cases. 



3. Direct force 

In the direct force case the situation is more complex 
than in the curvature force case and depend on the sign 
of both (Tth and (Ta. We will also distinguish between 
the long and short active correlation time as previously 
introduced. 



a. Long correlation time 

Examples of the fluctuation spectrum for the different 
limits of the tension have been plotted in (a) , (b) and (c) 
of Fig. [6] In this case for ath S> 1 (and hence (Ta > 0) 
the last term in equation ( S2 ) is always decreasing in I 



and we conclude that for ath S> 1 the active fluctuations 
are larger than the passive ones at low I and smaller at 
larger I. We stress here that this increase of the fluctua- 
tions at small I is present even though the active tension 
increases. Again we can see this by looking at the I — !■ oo 
limit 



1 



1 + Xa 



1 



(S13) 



The part proportional to Xa is quickly damped as I in- 
creases and the fluctuation spectrum goes under the pas- 
sive one as (Ta > ath whereas for small I the active noise 
directly increases the fluctuations. 
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In the case ath ~6 the small Xa expansion tells us 
that the active fluctuations are larger than the passive 
ones at Z = 2 and smaller for all larger I. In (a) of Fig. 
[6] we have put the fluctuation spectrum for a negative 
(Tth- In order to interpret this graph we flrst note that 
when (Tth -> — 6 more excess area are stored in the I ~ 2 
mode. Switching on the activity will then reduce the 
excess area stored in the I = 2 mode and increase the 
area stored in larger I. One can say that the activity 
drag the fluctuations from the / = 2 mode toward larger 
I mode. 



b. Short correlation time 

Examples of the fluctuation spectrum for the different 
limits of the tension have been plotted in (d) , (e) and (f ) 
of Fig. |6] In this case the only case where the term in 
the parenthesis of (S4) is monotonous is the case where 



CTth < 0. Then ) always goes from negative to 

positive as I increases. We observe again that the activity 
drags the excess area toward larger I. 

In the case ath > the small Xa expansion gives that 
for the case 1 ^ ath ^ ifn there are two crossover values, 
both are given in Table [ll] With the parameters we used 
the fluctuation spectrum keep the two crossovers at large 
Xa albeit with an increased / value for the crossover, see 
(e). 

For th e case 1 ^ ^ Cth the second positive root 
of (SlOl becomes larger than and there is only one 
crossover. The active fluctuations are larger than the 
passive ones for I < 1^ and smaller for I > If.. The ex- 



istence second root of (SIO) is the reason we could not 



conclude on the behaviors of the fluctuations using the 
derivative of the term in the parenthesis of (S4). Then 



the function in the parenthesis of ) is not guar- 

anteed to be monotonous in the whole I range, see (f). 



Permeation force 



Long correlation time 



Examples of the fluctuation spectrum for the different 
limits of the tension have been plotted in (a), (b) and 
(c) of Fig. [T] In this case no general conclusion can be 
derived from the derivative of the term in the parenthesis 



cases. All cases are found to have one crossover with 
active fluctuations larger than the passive ones for / < Ic 
for small Xa, but they show different behaviors at large 

Xa- 

For the case ath ~^ —Q and small Xa niost of the fluc- 
tuations are stored in the I = 2 mode and we observe 
a similar behavior as in the curvature and direct force 
cases, i.e., for a very large excess area the activity seems 
to drag the excess area toward large I values, see (a). 

For the case 1 ^ txth <C we can see that a new 
crossover appears (the root of 5{\uim\^) become smaller 
than ) when Xa. increases to ~ 1 and the flrst crossover 
disappears (the root of ) become smaller than 2) 

when Xa reach ~ 10^ resulting in a very different behavior 
depending on the value of the amplitude of the noise Xaj 
see (b). We note that the small Xa behavior is similar 
to the direct force case as the activity directly increases 
the low wavenumber fluctuations, whereas the large Xa 
behavior is more similar to the curvature force case where 
the activity increases the large wavenumber fluctuations. 

For the case 1 ^ ^ Cth we see again a second 
crossover appearing as Xa reach 10^ but the flrst crossover 
do not disappear. Instead we get a large increase of the 
I = 2 amplitude as Xa increases further, see (c). Note 
that we can make the second crossover disappear by de- 
creasing the excess area (data not shown). 



b. Short correlation time 

Examples of the fluctuation spectrum for the different 
limits of the tension have been plotted in (d), (e) and 
(f) of Fig. [7| In this case the second term in (|S4|) is 



always increasing with I and we can conclude in the same 
way as for the curvature case. The active fluctuations 
will be lower than the passive fluctuations for low / < Ic 
and higher for high I > Ic- This case is similar to the 
short correlation time limit of the curvature case for the 
following reason: if we calculate the I — >■ oo limit of the 
active fluctuation spectrum, Eq. (16) in the main text, 
for the permeation force case we get 



1 



(l+XaT40. 



(S14) 



of (S3 1 so we resort to the expansion in Xa for all the 



The large wavelength fluctuations are again increased by 
the activity and the small I fluctuations are correspond- 
ingly decreased through the tension to compensate. 
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(e) ' if) ' 

FIG. 5: The fluctuation spectrum as a function of I witli logaritlimic scale on both axis for the curvature force case. The passive 
fluctuations are in blue and the active fluctuations are in red. For all the subfigures we used /„! = 100 and k — 10. Figures (a), 
(b) and (c) are for 3> tm' while figures (d), (e) and (f) are for t'^ . Figures (a) and (d) are for ath — >■ —6. Here we used 

A = 1 which gives ath — —5.6, Xa — 10 for (a) and Xa — 0.3 10"'^ for (d). Figures (b) and (e) are for 1 <^ ath ^ l^- Here we 
used A = 0.1 which gives ath — 1592, Xa — 10 for (b) and Xa ~ 1.5 10~^ for (e). Figures (c) and (f) are for 1 <^ <^ ath- 
Here we used A = 10""^ which gives ath — 5 10^, Xa — 10'' for (c) and Xa — 1-5 10~^ for (f). Inset of the figure, is a density 
plot of the S{\uim\^) with the In/ along the horizontal axis and the base 10 logarithm of the amplitude of the activity in vertical 
axis. The graph is red when the active fluctuations are larger than the passive ones and is blue when the passive fluctuations 
are larger. The more light the color the larger is the absolute value of S{\uim\'^}- 
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Cej ' (f) ' 



FIG. 6: The correlation as a function of Z in a log log plot for the direct force case. We used the same color coding and 
parameters as in Fig. 1 except for: (a) Xa ^ 10^ (b) Xa 10^ (c) Xa ==; 10^ (d) Xa 10^ (e) Xa 10^ and (f) Xa 10^. 
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Cej ' (f) ' 

FIG. 7: The correlation as a function of Z in a log log plot for the permeation force case. We used the same color coding and 
parameters as in Fig. 1 except for: (a) Xa ==; 10^ (b) Xa 10^ (c) Xa ^ 10^ (d) Xa ==; 10-^ (e) Xa ^ 10"^ and (f) Xa 10"^ 



